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ABSTRACT 

We present two-dimensional Ha velocity fields for 20 late- type, disk-dominated spiral 
galaxies, the largest sample to date with high-resolution Ha velocity fields for bulgeless 
disks. From these data we derive rotation curves and the location of the kinematic cen- 
ters. The galaxy sample was selected to contain nucleated and non-nucleated galaxies 
(as determined from prior HST imaging), which allows us to investigate what impact 
the gas kinematics in the host disk have on the presence (or absence) of a nuclear star 
cluster. In general, we find that the velocity fields span a broad range of morpholo- 
gies. While some galaxies show regular rotation, most have some degree of irregular 
gas motions, which in nearly all cases can be either attributed to the presence of a 
bar or is connected to a rather patchy distribution of the Ha emission and the stellar 
light. There appears to be no systematic difference in the kinematics of nucleated and 
non-nucleated disks. Due to the large fields of view of the integral field units we use, 
we are able to observe the flattening of the rotation curve in almost all of our sample 
galaxies. This makes modeling of the velocity fields relatively straight- forward. 

Due to the complexities of the velocity fields, we obtain reliable determinations of 
the kinematic center for only 6 of our 20 sample galaxies. For all of these the locations 
of the nuclear star cluster /photometric center and the kinematic center agree within 
the uncertainties. These locations also agree for 7 more objects, despite considerably 
larger uncertainties as to the accuracy of the kinematic center. If we disregard all 
kinematically irregular galaxies, our study concludes that nuclear star clusters truly 
occupy the nuclei, or dynamical centers, of their hosts. Our results are thus consistent 
with in-situ formation of nuclear star clusters. Yet, many well- motivated formation 
scenarios for nuclear clusters invoke off-center cluster formation and subsequent sink- 
ing of clusters due to dynamical friction. In that case, our results imply that dynamical 
friction in the centers of bulgeless galaxies must be very effective in driving massive 
clusters to the kinematic center. 

Key words: galaxies: bulges; galaxies: nuclei; galaxies: star clusters; galaxies: spiral; 
galaxies: kinematics and dynamics 



1 INTRODUCTION 

What defines the center of a galaxy? This question is not 
merely academic, because throughout the last decade, a 
number of studies have found tight correlations between the 
global pro perties of g alaxies and the proper ties of their nu- 
clei (e.g. jFerrarese Merritt .2000.. : Gebhardt et al.1 |200QI : 
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iMarconi Huntil2003l : iHaring Ril3l2004l ). These global- 
to-nuclear scaling relations can be interpreted such that the 
mass assembly of a galaxy is pre-determined by its nuclear 
properties, or alternatively, that galaxy nuclei evolve in a 
way that is governed by the assembly of the entire galaxy. 

Either way, characterizing the nuclear properties has 
become an important diagnostic tool in constraining the for- 
mation mechanism(s) of galaxies. The question of where the 
galaxy nucleus - i.e. its center of mass - is located seems ob- 
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vious in ellipticals and bulge- dominated spirals, where very 
often a luminous active galactic nucleus (AGN) marks the 
location of a super-massive black hole (SMBH) which almost 
certainly marks the bottom of the potential well. 

However, in the latest Hubble-types, i.e. in bul- 
;eless, "pure" disk galaxies, (luminous) AGN are rare 
Satvapal et al.l [2009). and it is less obvious whether the 
galaxy disk rotates around a nucleus that follows in any 
way the scaling relations mentioned above. A number of re- 
cent studies have suggested that in late-type s pirals, the nu- 



cleus is marked by a massive stellar cluster ( Phillip s et al 
" Carollo et aLlflQQS : M atthews et al. 1999: Boker et al , 
20041 ). Such Nuclear Clusters (N Cs) are also present 
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in earlier type galax ies (e.^. rece ntly iBalcells et al.l l2003l : 
iLotz et al.i r2004: Cot e et al.1 12OO6I ). but the exact relation 
between NCs in galaxies of different Hubble types remains 
unclear to date. NCs have masses of - 10^ - 10'^ Mo 
(jWalcher et al. l2005l) and show stellar populations of multi- 
ple ages ("Walcher et al.l l2006l : iRossa et al.l [20061 : ISeth et al.1 
[2006 ), pointing towards them having a complex formation 
history. 

On average, the location of NCs appears to coincide 
with th e photometric cen ter (PC) as derived from isopho- 
tal fits ( Boker et al.l l2002l. However, the often irregular and 
asymmetric shape of late-type disk galaxies causes rather 
large uncertainties in defining the PC, and doubts have been 
raised on whe ther NCs actually define the bottom of the 
poten tial well ([Matthews Gallagherlliool : I Andersen et al.l 
I2OO8I ). 

Settling this question is important in order to rule 
out a number of suggested formation mechanisms for NCs. 
For example, if migration and/or merging of massive clus- 
ters is the dominant for mation mec hanism of NCs, as sug - 
gest ed bvlBekkil (l2007l)JCaDuzzo-Dolcetta fc^iocchil (|2008l ) 



and lAgarwal Milosavlievid (|2010l l. one would expect to 
find a number of NC s displaced from the nucleus, as e.g. 
iGeorgiev et~all (|2009l ) find for dwarf irregular galaxies. 

Since not all late- type spirals harbor an obvious NC 
([Boker et al.1 12OO2I ). one may also ask whether there is a 
galaxy property that prevents the formation of a NC. Nu- 
clear star formation may be suppressed in galaxies with ir- 
regular gas kinematics, and hence there may be systematic 
differences in the gas rotation patterns of galaxies with and 
without NCs. 

In order to address both these questions, we have ob- 
tained two-dimensional velocity fields of the ionized gas (as 
traced by the Ha line) for a sample of 20 late-type spiral 
galaxies. The use of Ha as a tracer for the general gas kine- 
matics is required for an accurate comparison of the posi- 
tion of the kinematic center (KG) and the PC, because it 
can be observed at high spatial resolution. While Ha may 
be affected by stellar winds, supernovae, or other deviations 
from the pure disk rotation, in general, it has been found to 
represen t well the overall ro tation field of the neutral (i.e. 
HI) gas (|S waters et al.ll2009l ). 

Similar studies have been performe d in the past (e.g. 
iGanda et al. ] l2006l : iBershadv et al.1 1 20 id) but this study is 
the first to focus exclusively on bulgeless spirals. In addi- 
tion, our analysis includes parameterised descriptions of the 
Ha rotation curve for all galaxies, thus enabling a direct 
comparison with the kinematics of other gas components. 

This paper is structured as follows: following this intro- 



duction, we describe in § [2] the galaxy sample, the obser- 
vations, and data reduction methods. In § [3l we detail the 
methods used to fit the Ha velocity fields, to derive the Ha 
rotation curves, and to extract the location of the KG. We 
discuss our results in § 2] and conclude in § [5] 



2 GALAXY SAMPLE AND OBSERVATIONS 

2.1 Sample Selection 

The galax y sample di s cusse d here was selected from the 
sample of iBoker et al.' (2002), and thus consists of spirals 
with late Hubble-type (Scd or later) and low inclination 
(< 40 ). In order to gauge the importance of a NC for 
the Ha kinematics, we selected galaxies with and without 
a NC in roughly equal parts. Given the visibility constraints 
from Calar AltcQ and WIYPsQ, we identified a sample of 20 
objects, summarized in Table [1] The object distances are 
between 6 and 28 Mpc, which implies that the field of view 
of the Ppak and SparsePak instruments cover between 2 
and 9.5 kpc in radius. The main instrumental parameters of 
Ppak and SparsePak are summarised in Tabled 

2.2 Ppak Data 

During the nights of May 7th and 8th, 2007, thirteen galax- 
ies of our sample were observed at the 3.5 m telescope of 
the Calar Alto observatory, using the Potsdam Multi Aper- 
ture Spectrograp h (PM AS, Roth et al. 2005) in its Ppak 
mode (|Kelz et al.ll2006l ). The Ppak science fibre bundle con- 
sists of 382 fibres of 2.7 arcsec diameter each, of which 331 
(the science fibres) are concentrated in a hexagonal bundle 
covering a field-of-view of 72" x 64" with a filling factor of 
^^65%. The sky background is sampled by 36 additional fi- 
bres, distributed in 6 bundles of 6 fibres each, distributed 
along a circle ^ 90 arcsec from the center of the instrument 
FOV. The sky- fibres are distributed among the science fi- 
bres within the pseudo-slit in order to have a good charac- 
terization of the sky; the remaining 15 fibres are used for 
calibration purposes. Cross-talk between adjacent fibres is 
estimated to be less than 5% when using a simple aperture 
extraction (Sanchez 2006). Adjacent fibres in the pseudo-slit 
may cover very different locations on the sky, thus further 
reducing the effect of cross-talk. 

The J1200 grating, mounted backwards in 2nd order, 
was used for all observations. It covers the wavelength range 
-6350-6690 A with a spectral resolution of FWHM -0.5 
A(R —11000). For each galaxy, two exposures were taken, 
with exposure times between 900s and 1500s, depending on 
the target brightness. The nights were clear, with a slightly 
elevated extinction (Ay —0.18 mag), but stable in both 
cases. The seeing was variable, ranging between 0.8^^ and 
1.3^^ A spectrophotometric standard star was observed each 
night, in order to correct for the transmission curve of the 



^ The German-Spanish Astronomical Center, Calar Alto, is 
jointly operated by the Max- Planck- Inst it ut fr Astronomic Hei- 
delberg and the Instituto de Astrofsica de Andaluca (CSIC). 
^ The WIYN Observatory is a joint facility of the University of 
Wisconsin-Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatories. 
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Table 1. Galaxy Sample 



Galaxy 


Hubble type 


RA 


DEC 


D 


Li 






Mo 






(lih:mm:ss.ss) 


(dd:mm:ss.s) 


(Mpc) 


(109 Lo) 


(km s-l) 


(mag) 


(mag/arcsec' 


UGC 35742 


SAcd 


06:53:10.4 


+57:10:40 


20.7+3.0 


10.0 


167 


-11.90+0.18 


18.9+0.1 


NGC25522 


SAm 


08:19:20.5 


+50:00:35 


11.9+2.0 


2.3 


62 


-12.04+0.01 


20.15+0.1 


UGC 44992 


SABdm 


08:37:41.5 


+51:39:09 


12.8+2.0 


0.23 


51 


-8.59+0.63 


19.85+0.05 


UGC 5288^ 


Sdm 


09:51:17.0 


+07:49:39 


6.0+1.2 


0.072 


49 


- 


19.8 ±0.1 


NGC 32062 


SBcd 


10:21:47.6 


+56:55:50 


20.7+0.4 


3.0 


78 


- 


18.78±0.05 


NGC33461 


SBcd 


10:43:38.9 


+14:52:19 


22.4+3.0 


8.8 


124 


-11.78 ±0.01 


18.3±0.2 


NGC 34232 


SAcd 


10:51:14.3 


+05:50:24 


11.3+0.6 


1.4 


127 


-11.84+0.05 


17.9±0.15 


NGC3445I 


SABm 


10:54:35.5 


+56:59:26 


17.5+2.0 


3.4 


148 


-13.42+0.10 


18.1±0.1 


NGC4204I 


SBdm 


12:15:14.3 


+20:39:32 


7.9+2.0 


0.25 


- 


-10.26+0.02 


19.6±0.1 


NGC4299I 


SABdm 


12:21:40.9 


+11:30:12 


16.8+2.0 


3.9 


109 


-11.73+0.04 


19.0±0.5 


NGC 4496ai 


SBm 


12:31:39.2 


+03:56:22 


15.6+1.2 


5.8 


94 


-11.99+0.02 


18.9±0.1 


NGC4517ai 


SBdm 


12:32:28.1 


+00:23:23 


24.5+3.0 


1.4 


69 


- 


19.8±0.1 


NGC 45402 


SABcd 


12:34:50.8 


+15:33:05 


16.8+2.0 


4.8 


83 


-12.29+0.02 


18.05±0.05 


NGC46251 


SABm 


12:41:52.7 


+41:16:26 


8.2+2.0 


0.51 


39 


-10.61+0.08 


17.5±0.2 


NGC4904i'2 


SBcd 


13:00:58.6 


-00:01:40 


20.0+2.0 


2.1 


105 




17.3±0.2 


UGC85161 


Scd 


13:31:52.6 


+20:00:04 


20.4+2.0 


1.6 


60 


-10.97+0.09 


18.7±0.1 


NGC 56692 


SABcd 


14:32:43.5 


+09:53:26 


18.0+2.0 


8.6 


98 


-10.03+0.01 


18.5±0.1 


NGC 5789i'2 


Sdm 


14:56:35.5 


+30:14:03 


33.0+3.0 


4.3 


123 




19.9±0.05 


NGC 5964^ 


SBd 


15:37:36.3 


+05:58:24 


26.5+2.0 


21.6 


121 


-12.62+0.06 


18.9±0.1 


NGC6509I 


Sd 


17:59:25.3 


+06:17:13 


28.2+3.0 


6.3 


218 


-13.08+0.07 


17.8±0.1 



Note. - Ga laxies with ^ were observe d with CAHA/Ppak, those with 2 with WIYN/SparsePak. RA, DEC and Hubble types are from 
the RC3 (jde Vau couleurs et al. 1991) through NED. Distances have been assembled through NEDs summary statistics and are from 
redshift independent measurements, where possible. If only one measurement was available we assumed a default uncertainty of 
max(2Mpc, D*0.1). For UGC 8516 and NGC 5789 the distance co mes from the redshift, corrected for Virgo, Great Attractor and 
Shapley cluster infall, following the formulae of iMould et al ] (|2000l V Li and Vc have been assembled through Leda and use their itc 
and Vmax param eters. Central surface brightnesses /j,o have been derived from the HST I-band surface brightness profiles (SBP) of 
iBoker et alJ (|20Q2h . For galaxies without a NC the SBPs are flat in all cases and /j^q is well defined. For galaxies with a NC it is unclear 
to what extend the disk extends further into the NC, we therefore chose to use the surface brightness just outside the NC as a measure 
of /io- Error bars give our estimate of the systematic uncertainty of this number. 



instrument. Note though that spectrophotometric accuracy 
is not required for our analysis. 

The data were reduced using R3D ( Sanched l2006r). 
in combination with irafEI packages and E3D (|Sanchezl 
The reduction consists of the standard steps for fibre- 
based integral-field spectroscopy. A master bias frame (cre- 
ated by averaging all bias frames observed during the night) 
was subtracted from the science frames. Exposures of a given 
sky position were median-combined using IRAF routines, 
thus clipping any cosmic rays. The locations of the spectra 
on the CCD were traced using an exposure of a continuum 
lamp taken before the science exposures. Each spectrum was 
then extracted from th e science frame s, and stored in a row- 
stacked-spectrum file (|Sanchej[200i l. Wavelength calibra- 
tion was performed using the position of Ne lines in lamp 
exposures obtained before and after each pointing, yielding 
an accuracy of rms~ 0.15 A. Differences in the relative fibre- 
to-fibre transmission throughput were corrected by compar- 
ing the wavelength-calibrated science frames with the cor- 
responding frames derived from sky exposures taken dur- 
ing twilight. Then, the data was corrected for the average 
instrument throughput curve, by comparing the observed 
spectrum of the spectrophotometric calibration star with 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



Table 2. Instrumental parameters 



Ppak SparsePak 



# science fibres 


331 


75 


# sky flbres 


36 


7 


size of flbres 


2'.'7 


4('675 


fleld-of-view 


72'' X 64'' 


70" X 70" 


fllling factor 


-57% 


~ 27% 


wavelength range 


6350 - 6690A 


6500 - 6900A 


instrumental FWHM 


O.5A 


0.65A 



a flux-calibrated one. Finally, a contemporaneous (average) 
night sky spectrum was obtained by combining the spectra 
of the 36 sky fibers, and subtracted from the science spectra. 
The relative location of the final science spectra on the sky 
were obtained via the standard Ppak position table. 

2.3 SparsePak Data 

Nine galaxies in our sample were observed using SparsePak 
(Bershady et al. 2004) on the 3.5m WIYN telescope during 
the course of four nights, namely March 27-30, 2007. 

SparsePak is a fiber optic-array that feeds light from 
the WIYN Nasmyth f/6.3 focus imaging port to the Bench 
Spectrograph. The SparsePak IFU consists of 82 fibers with 
a diameter of 4.675^^ each, arranged in a sparse grid with a 
field of regard of ^ 70^^ x 70^^ thus offering good coverage 
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and sampling of our target galaxies. The spectrograph was 
configured with the Bench Spectrograph Camera (BSC) and 
an Echelle grating with 316 lines/mm, used in 8th order. It 
covers the wavelength range 6500 A < A < 6900 Awith a dis- 
persion of 0.195 A/pix (8.6 km s~^/pix) and an instrumental 
FWHM of 0.67 A (30.5 km s"^). 

This rather high resolution is mandatory for our pur- 
pose: fitting velocity field models to galaxies with observed 
rotation velocities of roughly 100 km s~^ requires a cen- 
troiding accuracy of 5 km s-^The BSC images the spec- 
trograph onto a T2KA CCD with 2048x2048 pixels. The 
spectra are aligned along the columns of the CCD. The chip 
has a read noise of 4 e~ and was used with the standard 
gain of 2.1 e~/ ADU. The system throughput for this setup 
is roughly 4% (jBershadv e t al. 2005). 

The sparse grid can be filled in with 3 pointings, and 
we observed 8 of the 9 galaxies with 3 pointings (NGC 5789 
was only observed with two pointings). When WIYN first 
points to a target, the slit- viewing camera is used to put 
the PC roughly coincident with fiber 52 at the center of 
the SparsePak array (if the target has a surface brightness 
that is too low, as is the case for many of these targets, we 
trust that the WIYN pointing is accurate enough to deliver 
the science target onto the fiber array after offsetting from 
a nearby star). After 2 x 20 minute exposures (we use 2 
exposures to be able to better reject cosmic rays), guiding 
is stopped, and the telescope is offset by 5.6^^ towards the 
South and guiding is resumed. After the next 2 x 20 minute 
exposures, guiding is paused, the telescope is offset by 4.9^^ 
West and 2.8^^ North from the second position. Guiding is 
resumed, and the final 2 x 20 minute exposures for a galaxy 
are taken. 

Data were overscan- and bias-corrected and trimmed 
using the NOAO IRAF package ccdproc. Cosmic ray re- 
jection was performed b efore spectral e x tracti on, using 
the method described in lAndersen et al.l (|2006l V Follow- 
ing cosmic-ray cleaning, basic spectral extraction, flatten- 
ing, and wavelength calibration were done using the IRAF 
package dohydra. During this process, we made use of bias 
frames, dome flats, and Thorium Argon emission spectra 
that were taken each night. Finally, the emission from the 
night sky was subtracted by averaging the spectra of the 
seven sky fibers, and subtracting the result from each of the 
75 source spectra. 



2.4 Image Registration 

In order to correlate the results of the kinematical analysis 
with the galaxy morphology, and in particular the positions 
of KC, PC, and NC, we use archival HST/WFPC2 F814W 
images (roughly corresponding to the Joh nson I-band) from 
the snapshot survey of Boker et al.l For each galaxy, 

we first created a continuum image from the IFU data cube 
by fitting the continuum in a spectral window free of emis- 
sion lines. This should allow a fair comparison to the I-band 
images. We then smoothed the HST image with a Gaussian 
beam whose width matches the seeing-limited resolution of 
the IFU observations (typically 1.0^^). 

We cross- correlated both images as follows: for a given 
position of the (smoothed) HST image, we extracted its flux 
within the footprint of each IFU spaxel. We then fit a linear 
relation between the IFU and the HST spaxel fluxes and 
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Figure 2. Ho; emission lines from two Ppak spaxels around the 
center of NGC6509 comparing the double-Gaussian fits (top) to 
the single-Gaussian fits (middle). The bottom panels show the fit 
residuals to the fit in the middle panel. Ho; was detected in 3078 
of 4303 Ppak spaxels and 1220 of 1950 SparsePak spaxels. Most 
lines can be well fit by a single Gaussian, as shown in the plots. 



tabulated x^- We repeated this process over a grid of offsets 
between both images, thus "mapping" • The offset grid 
has a granularity corresponding to five times the pixel size 
of the WFPC2 data, i.e. 0.232''. The location of the mini- 
mum within the offset grid was used as the best registration, 
and the minimum value to estimate the uncertainty in 
the registration (see Fig. [T] for an illustration). 



2.5 Line Fitting and Extraction of Velocity Fields 

Once spectra were processed as described in § [2l we identi- 
fied Ha emission-lines and measured fluxes, widths, centers 
and the corresponding errors for lines in a given spectral 
window. We detected Ha emission lines and measured their 
fluxes, widths and centroids in 3078 of 4303 Ppak and 1220 
of 1950 SparsePak object spectra (72% and 63% detection 
rate, respectively). While the galaxy described in Andersen 
et al. (2008), NGC 2139, contained multiple kinematic com- 
ponents in each spectra, most emission lines in this sample 
were best fit by a single Gaussian line. Furthermore, typi- 
cally 10 spectra per galaxy were best fit by two Gaussians 
with centers that were close to being coincident, i.e., the 
lines appeared to have a strong core with broader, low-level 
wings. Still, the center positions and widths of these lines 
can be very well fit by a single Gaussian (see Fig. [2])- Only 
about 2-3 lines per galaxy exhibited two de-coupled Gaus- 
sian profiles. Almost all of these double-line features occur 
at the very center of the galaxies. In case of double line fea- 
tures, we used the more highly peaked, dominant component 
of the line fit to derive the velocity map. 
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Figure 1. Registration of the IFU data to the HST images: Left: HST F814W image convolved with the SparsePak PSF and sampled 
to the SparsePak pixelsize in comparison with the continuum levels from the SparsePak spaxels (second panel) for NGC3423. Third 
panel: Spaxel-ensquared HST WFPC2 F814W flux vs. SparsePak continuum flux with best-flt linear regression through (0,0). Right: 
Error contours of the image IFU registration based on the map, indicating an astrometric accuracy of the SparsePak data of ±0(^3. 



3 ANALYSIS 

3.1 Kinematic Modelling 

In order to interpret the Ha velocity fields that result from 
the analysis described in the last section, we fit them with 
a kinematic model. While there are arguments for fitting a 
more flexible tilted ring model, which allows for warps in 
the gas disks, the sometimes sparse sampling of the over- 
all velocity field of our data did not allow this approach 
for all galaxies. In the interest of a uniform analysis for 
th e entire sample, we therefore chose to follow the approach 
of lCourteaul (|l997l l who proposed the following parametric 
form for the (de-projected) rotation curve: 

vir) = + — i;carctan(^^^ — ^) (1) 
TT rt 

ICourteaul (|l997l l shows that this function describes the 
rotation curves of spiral galaxies well. Its specific advantage 
in our context is that it has only two intrinsic free param- 
eters, i.e. scale radius and velocity scale, which is the min- 
imum number of parameters possible for a rotation curve. 
The modeled two-dimensional velocity fields thus have seven 
free parameters: the systemic velocity VQ^ihe asymptotic ve- 
locity Vc^ the center positions xo and ^/o, the scale radius rt, 
the inclination i and the position angle PA. 

The fits are carried out via a custom-written Markov 
Chain Monte Carlo code. In addition to the best fit parame- 
ters, w e also obtain their uncertainties. As in Andersen et ahl 
(|2008h we find that the formal, statistical errors on the Ha 
velocities are too small to yield a useful measure of the qual- 
ity of the fit. We thus add a second, additive error term a mod 
to make sure that the reduced of the best fit is not un- 
reasonably high. We found that although our best-fit center 
does not depend on (Jmod, the uncertainties in the position 
of the KC do. Also, this ad-hoc additional error term, ham- 
pers our capability to independently assess the quality of 
the final fit. We chose to set cTmod = 8 kra/s for all galax- 
ies (compare e.g. iKamphuis k, Sancisill 19931 : 1 Andersen et al.l 
l2006l : ISellwood Sanchezll201Ql . and references therein). 

The best-fit kinematic models are presented in column 4 
of Fig. m and the derived rotation curves, including the best 
fit to Eq.[T] are shown in Figure[5] The Figures show that our 




Ledo: opporent moximum rotation velocity 



Figure 3. The projected maximal rotation velocity of HI gas 
from the Hyper leda database against the projected velocity of the 
model at 3 times the scale radius rt- We find a very acceptable 
correlation. The line indicates the expected one-to-one correla- 
tion. 



model velocity fields seem to be a fair representation of the 
observations, confirming that also the latest Hubble types 
can be fit reasonably well with this simple functional form. 
The resulting best-fit kinematic parameters are summarised 
in Table El 

Since we are mostly interested in the position of the KC, 
we have checked the dependence of the fit results on the cho- 
sen starting values for the center position, for i, PA, vq^ etc. 
We find that the fit always converges to the same solution, 
within ^30% of the nominal error bar, demonstrating that 
our results are insensitive to the choice of starting values. 
The derived coordinates of the KCs are given in Table [H 
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Table 3. Final Kinematic Parameters 



Galaxy 




Vsys,NED 


PA 




i 




Center RA 


Center DbC 


Center error 


^circ 








/I —1\ 

(km s j 


/I —1\ 
(km s j 


( ) 




^o\ 
( ) 




(nn:mm:ss.ssj 


(dd:mm:ss.sj 


(arcsecj 


/I — 1 

(km s 


) 


arcsec 


(1) 


(2) 


(3) 


/ /I \ 

(4) 




(5) 




(6) 


(V 


(8) 


(9) 




(10) 


UGC 35742 


1438±2 


1441 


98.4±0, 


,5 


25.3±0, 


,2 


06:53:10.43 


+57:10:37.9 


0.9 


167.2+0, 


.3 


198.4+0.4 


NGC 25522 


511±3 


524 


62.4±0, 


,9 


25.4±0, 


,4 


08:19:19.44 


+50:00:48.7 


2.3 


62.2+1, 


,6 


221.0+1.9 


UGC 44992 


685±2 


691 


141. 9±1, 


.7 


29.0±0, 


.2 


08:37:41.49 


+51:39:13.6 


6.8 


50.9+1, 


.6 


200.2+2.7 


UGC 5288^ 


576±6 


556 


56.7±0. 


,5 


29.9±0. 


,2 


09:51:16.88 


+07:49:47.3 


4.3 


48.6+0, 


.4 


179.0+0.4 


NGC 32062 


1151±1 


1150 


183. 0±0, 


,5 


26.5±0, 


.1 


10:21:47.81 


+56:55:49.8 


0.6 


78.4+0, 


.4 


229.3+1.2 


NGC33461 


1274±1 


1260 


-67.1±0, 


,6 


29.4±0, 


,3 


10:43:38.85 


+14:52:16.6 


0.6 


123.6+0, 


,8 


207.1+1.3 


NGC 34232 


1004±2 


1011 


45.1±0, 


,8 


19.0±0, 


,2 


10:51:14.36 


+05:50:24.3 


1.1 


127.0+0, 


.5 


295.6+1.2 


NGC3445I 


2048±1 


2069 


-9.8±1, 


,1 


27.6±0, 


,4 


10:54:35.21 


+56:59:23.7 


1.9 


148.1+2, 


,3 


132.1+4.9 


NGC4204I 


870±3 


856 


240.2±2, 


,5 


13.6±1, 


,7 


12:15:14.36 


+20:39:29.6 


2.6 


50.0+2, 


.8 


238.0+19.3 


NGC4299I 


237±2 


232 


-110.3±1, 


,0 


28.7±1, 


,2 


12:21:40.59 


+11:30:11.4 


0.9 


109.1+1, 


,1 


189.9+4.8 


NGC 4496a^ 


1747±3 


1730 


47.5±0, 


,3 


30.4±0, 


,1 


12:31:39.75 


+03:56:18.9 


1.7 


94.3+0, 


.2 


221.3+0.5 


NGC4517ai 


1525±2 


1509 


-156. 6±1, 


,5 


33.7±1, 


,2 


12:32:28.32 


+00:23:28.7 


1.9 


68.6+0, 


,8 


175.5+5.8 


NGC 45402 


1291±2 


1286 


13.9±1, 


,7 


27.9±1, 


,1 


12:34:50.91 


+15:33:06.3 


1.4 


83.4+1, 


,5 


201.0+2.6 


NGC46251 


621±1 


609 


-55.8±1, 


,1 


13.3±0, 


,3 


12:41:53.06 


+41:16:24.0 


1.3 


38.7+0, 


.2 


188.3+1.3 


NGC4904I 


1180±1 


1189 


-134.0±0, 


,5 


38.5±0, 


,6 


13:00:58.62 


-00:01:37.8 


0.3 


105.2+0, 


.2 


235.1+1.8 


NGC 49042 


1162±2 


1189 


227.2±0, 


,9 


39.0±0, 


,5 


13:00:58.54 


-00:01:37.9 


0.4 


105.2+0, 


.9 


247.5+2.8 


UGC85161 


1026±3 


1023 


14.5±0, 


,7 


43.3±0, 


,3 


13:31:52.60 


+20:00:03.9 


1.4 


60.2+0, 


.9 


154.5+1.2 


NGC 56692 


1368±2 


1371 


69.1±0, 


,7 


35.6±0, 


,0 


14:32:44.06 


+09:53:29.5 


0.5 


98.4+0, 


.5 


186.6+0.6 


NGC57891 


1811±1 


1805 


150.3±0, 


,8 


27.4±0, 


,1 


14:56:35.58 


+30:14:02.4 


1.2 


122.9+1, 


.0 


140.3+.6 


NGC 57892 


1809±2 


1805 


151. 7±2, 


,9 


27.0±0, 


,2 


14:56:35.56 


+30:14:01.1 


2.2 


122.9+1, 


,3 


148.8+4.5 


NGC 59641 


1457±4 


1447 


131. 8±0, 


,9 


38.4±0, 


,2 


15:37:36.94 


+05:58:17.3 


5.1 


120.8+0, 


.8 


220.1+1.2 


NGC6509I 


1780±1 


1813 


-81.6±0, 


,4 


49.3±0, 


,7 


17:59:25.46 


+06:17:10.5 


0.3 


218.1+0, 


,4 


213.7+0.7 



Note. - Galaxies with ^ were observed with CAHA/Ppak, those with 2 with WIYN/SparsePak. The kinematic parameters as derived 
in this paper are (2) systemic velocity , (4) Position to the major axis, (5) inclination, (6) and (7) RA and DEC of the KG, (8) error on 
the kinematic center, (9) Circular velocity, defined as the maximum velocity in the rotation curve fits, and (10) the scale radius of the 
rotation curve. For comparison to the modelled systemic velocity (2), we give the systemic velocity assembled through NED (3). 



3.2 Uncertainties in the kinematic modelling 

For an independent check, we compared our results to the 
apparent maximal rotation velocity of the gas (parame- 
ter v maxq) from the Hyp erLEDA d atabasj^ ([Paturel et alJ 
I2OO3I ) . As pointed out by ICourteaul (| 19971 ), the parameter 
Vc in Eq. [T] is in itself not a good measure of the maximal 
velocity. Instead, we use V3t, i-e. the model velocity at three 
times rt . Figure [3] shows that vst correlates well with the 
Vmax parameter from LED A, after correction for the galaxy 
inclination as derived from our fits. There are a number of 
uncertainties underlying the plotted values giving rise to the 
scatter in the plotted relation, like irregularities in the ve- 
locity fields for both HI and Ha or the fact that we do not 
sample the full velocity field with our Hadata. Given these 
uncertainties, the scatter is surprisingly low. This demon- 
strates the reliability of our kinematic model^. 

A further consistency check can be performed by 
comparing the dynamical models derived from Ppak and 
SparsePak data for NGC 4904 and NGC 5789, the two 
galaxies which have been observed with both instruments. 
As can be seen from Fig. the observations agree well in 
both cases. The fit results also agree well, both in the loca- 
tion of the KG, and in the overall shape of the contours, 
despite marked differences in the spatial sampling. We take 



this as confirmation of both the quality of our data and the 
robustness of our analysis method. 

However, our main uncertainty in the velocity field 
modelling is the actual complexity of the velocity fields. Our 
simple model does e.g. not include the effects of streaming 
motions due to bars or the effects on the velocity field due to 
interactions with a companion galaxy. We have performed 
an eye-ball check of our trust in the velocity field modelling, 
by looking at the spacing of the contours of the fit to 
the dynamical centre (see Fig. |4] column 5), weighing in the 
presence/absence of a bar /companion, the filling factor of 
Halpha emission, the overall regularity of the velocity field 
and the quality of the fit to the rotation curve. On that basis 
we have assigned quality flags from (unusable) over 1 (low 
trust) to 2 (well modelled) to each velocity field. This quality 
assessments number is added as a column in Table |4l This 
exercise has been somewhat sobering as to our ability to 
systematically determine dynamical galaxy centers, as only 
7 velocity fields out of 22 obtain the label "well modelled" . 
Only two of these are neither barred nor show any other sign 
of problems (NGC 3423, NGC 3206). Late type disk galaxies 
tend to be irregular and a well-defined center simply does 
not exist in quite a few of them. The implications of this 
will be further discussed in Section U 



^ http://leda.univ-lyonl.fr 

^ We note that it is doubtful whether further exploration of our 
data for Tully-Fisher relation purposes is reasonable, due to the 
many kinematic uncertainties we identify. 



3.3 Comparing kinematic centers and 
photocenters 

Following image registration as described in § 12.41 we can 
measure the absolute coordinates of the best fit KG. Ta- 
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Figure 4. Montage of the sample galaxies sorted by right ascension, one line per galaxy. Panels from left to right: (i) SDSS r-Band (when 
available) or DSS R-Band image, with two footprints per galaxy overplotted. The small square footprint indicates the HST WFPC2 
field-of-view, the larger square gives the SparsePak and the hexagon the Ppak footprint, (ii) Ha flux map and (iii) observed Ha velocity 
map in comparison to (iv) the modelled velocity field (masked with the observed fiux map). The right panel shows the HST WFPC2 
F814W image with a transparent overlay of the velocity map. Gap s in the v elocity map are filled in with model values. The white cross 
indicates the position and uncertainty of the PC (as derived by Bo ker et al.ll20 02V black contours indicate the position and uncertainty 
^^^iS^fiX'&^^SiA^if/^S §db j^^gj^^^^^^^ ^^^^ other parameters), derived in this paper. 
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Figure 4 - continued 
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Figure 4 - continued 
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Figure 4 - continued 
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Table 4. Compilation of Kinematic Centers, Nuclear Cluster Positions, and Photometric Centers 





Kinematic Center 


offset 


Nuclear Cluster Position 


offset 


Photometric Center 




Galaxy 


RA 


DEC 


(KC-NC) 


RA 


DEC 


(PC-KC) 


RA 


DEC 


(I 




hh:mm:ss.ss 


dd:mm:ss.s 




hh:mm:ss.ss 


dd:mm:ss.s 




hh:mm:ss.ss 


dd:mm:ss.s 




UGC 3574^ 


06:53:10.43 


+57:10:37.9 


1.7/170.6 


06:53:10.39 


+57:10:39.5 


5.7/ 572.0 


06:53:10.39 


+57:10:39.5 


1 


NGC25522 


08:19:19.44 


+50:00:48.7 


21.3/1228.9 


08:19:20.38 


+50:00:32.8 


20.6/1188.5 


08:19:20.28 


+50:00:32.4 


1 


UGC 4499^ 


08:37:41.49 


+51:39:13.6 


4.5/279.3 


08:37:41.43 


+51:39:09.2 


5.1/316.5 


08:37:41.44 


+51:39:08.6 


1 


UGC 5288-*- 


09:51:16.88 


+07:49:47.3 




- 


- 


8.8/256.0 


09:51:17.09 


+07:49:39.1 





NGC 3206^ 


10:21:47.81 


+56:55:49.8 




- 


- 


0.2/20.1 


10:21:47.79 


+56:55:49.7 


2 


NGC 3346^ 


10:43:38.85 


+14:52:16.6 


0.6/65.2 


10:43:38.84 


+14:52:17.2 


0.9/97.7 


10:43:38.82 


+14:52:17.4 


2 


NGC 3423^ 


10:51:14.36 


+05:50:24.3 


3 1 /1 69 8 


10:51:14.31 


+05:50:24.3 


3 1 /1 69 8 


10:51:14.31 


+05:50:24.3 


2 


NGC 3445^ 


10:54:35.21 


+56:59:23.7 


5 8/492 1 


10:54:35.55 


+56:59:26.4 


6 1 /517 5 


10:54:35.59 


+56:59:25.9 


1 


NGC 4204-*^ 


12:15:14.36 


+20:39:29.6 


1.4/53.6 


12:15:14.43 


+20:39:28.8 


1.4/53.6 


12:15:14.40 


+20:39:29.6 





NGC 4299^ 


12:21:40.59 


+11:30:11.4 


2.9 /236.2 


12:21:40.39 


+ 11:30:11.8 


2 3/187 3 


12:21:40.44 


+11:30:10.6 


2 


NGC 4496a-'^ 


12:31:39.75 


+03:56:18.9 


7.8/589.9 


12:31:39.24 


+03:56:20.5 


7.6/574.8 


12:31:39.26 


+03:56:20.9 


1 


NGC4517a^ 


12:32:28.32 


+00:23:28.7 




- 


- 


7.0/831.5 


12:32:28.15 


+00:23:22.2 


1 


NGC 4540^ 


12:34:50.91 


+15:33:06.3 


2.7/219.9 


12:34:50.89 


+15:33:06.9 


0.7/57.0 


12:34:50.89 


+15:33;06.9 


1 


NGC 4625^ 


12:41:53.06 


+41:16:24.0 


1.6/63.6 


12:41:53.06 


+41:16:25.6 


1.6/63.6 


12:41:53.06 


+41:16:25.6 


1 


NGC49041 


13:00:58.62 


-00:01:37.8 








2.8/271.5 


13:00:58.59 


-00:01:37.9 


2 


NGC 4904^ 


13:00:58.54 


-00:01:37.9 








0.8/77.6 


13:00:58.59 


-00:01:37.9 


2 


UGC85161 


13:31:52.60 


+20:00:03.9 


0.6/59.3 


13:31:52.56 


+20:00:03.8 


0.9/89.0 


13:31:52.54 


+20:00:03.6 


1 


NGC 56692 


14:32:44.06 


+09:53:29.5 


1.2/104.7 


14:32:44.09 


+09:53:30.6 


1.0/87.3 


14:32:44.09 


+09:53:30.4 


2 


NGC57891 


14:56:35.58 


+30:14:02.4 








4.0/640.0 


14:56:35.84 


+30:14:03.1 


1 


NGC 5789^ 


14:56:35.56 


+30:14:01.1 








4.8/767.9 


14:56:35.84 


+30:14:03.1 


1 


NGC 5964^ 


15:37:36.94 


+05:58:17.3 


15.0/1927.1 


15:37:36.11 


+05:58:25.9 


14.7/1888.6 


15:37:36.13 


+05:58:25.7 





NGC65091 


17:59:25.46 


+06:17:10.5 


2.5/341.8 


17:59:25.36 


+06:17:12.5 


2.6/ 355.5 


17:59:25.34 


+06:17:12.4 


1 



Note. - Galaxies with ^ were observed with CAHA/Ppak, those with ^ with WIYN/SparsePak. The position of the kinematic center 
as derived in this paper is given in columns RA^c and DECxc". The Nu c lear C luster coordinates (RAjvc/DECivc) and coordinates 

of the photometric center (RApc/DECpc) are taken from Bok er et al ] (l2002l V The offsets between nuclear cluster and kinematic 
centre (KC-NC) and photometric centre and kinematic centre (PC-KC) are given in arcseconds and parsecs (using the distances given 

in Table [T] The parameter g is a measure for how well the velocity fields can be modelled, where denotes a bad fit to the data, 1 

denotes models with low trust and 2 denote models of high fidelity. 



NGC5964 SDSS r-Band 




Figure 4 - continued 
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ble [4] summarizes the KC positions, along with those of 
the PC and NC (if present). The latter two were derived 
from isopho t al fits to the WFPC2 image, as described in 
iBoker et all (|2002l l. Ah coordinates in Table H refer to the 
J2000 coordinate system of the HST image, which has an 
absolute accuracy with respect to the ICRS system of ^ l" . 
In Fig. [6l we compare the relative positions of the NC and 
the KC as well as the PC and the KC for the non-nucleated 
galaxies. 

Generally, we find that there are only a few galaxies 
(6 out of 20) with a velocity field regular enough (trust 
level 2 in Table [4|) to meaningfully compare the locations 
of NC, PC and KC. For all of these galaxies these locations 
agree within the uncertainties, although some of the galax- 
ies show bars. On the other hand, for all but one galaxy 
with offset KC we can always identify a reason. NGC 2552, 
NGC4299, NGC 5789, have a very fuzzy or patchy appear- 
ance in Ha and in stellar emission. These galaxies may 
simply not be in regular rotation. UGC5288, NGC 4204, 
NGC 4496, NGC 5669, NGC 5964, NGC 6509 show a strong 
bar and the velocity field is clearly affected by large scale 
streaming motions along the bar. NGC 4625 is lopsided from 
a strong m=l mode, which seems to displace the KC. There 
are only two potentially odd cases: in NGC4517 there is a 
large offset between KC and NC. However, this galaxy prob- 
ably has a nuclear starburst (see Ha map) that affects the 
central part of the velocity field, although a possible bar is 
weak if existent. UGC 3574 looks very regular, yet NC and 
KC show an offset of the size of the 1 a error bar. However, 
in this galaxy the central ^" unfortunately do not show Ha 
emission, which might hamper our ability to determine an 
accurate KC. We are left with no confirmed offset between 
the NC and the KC. 



3.4 Rotation curves and velocity residuals 

Rotation curves of our sample galaxies are shown in Figure 
[5] As mentioned in Section [3.11 we find that the functional 
form given in Eq. [T] represents the rotation curves of latest 
type spirals well in general. To quantify this statement in 
terms of the scientific goal of our paper we have used the 
rotation curves of our best models (trust level 2 in Table [4|) 
and our bad models (trust level 0) . There is no obvious trend 
with nucleatedness (4 out of 6 well-fit galaxies are nucleated 
and 1 out of 3 not- well- fit galaxies). 

Those rotation curves which are least consistent with 
the arctan functional form are dominated by irregularities of 
their host galaxies. The irregular rotation curve of NGC 4625 
can be explained by an ongoing interaction, while NGC 5964 
is strongly barred, and UGC 4499 has a very sparsely sam- 
pled velocity field. 

The rotation curve shapes of l ate ty pe galaxies have also 
been discussed in I S waters et al.l (|2009l ) from HI data. The 
morphological variety of their rotation curves is similar to 



3.5 Dynamical Friction Timescales 

Based on Chandrasekhar's formula we calculate the dynam- 
ical friction timescales for star clusters in our galaxies using 



equation 8.12 of iBinnev &: Tremaind (|200 



tfri 



InA GM ' 



(2) 



where InK is the Coulomb logarithm, Vc the typical ve- 
locity of the stars in the galaxy at radius r^, and M 
is the mass of the test particle that is going through 
the much larger ma s s of t he underlying galaxy. Although 
Binnev & Tremain i (|2008') derived Eq. [2] for the decay of 
black hole orbits in a singular isothermal sphere, they explic- 
itly state that it is approximately correct even for mass dis- 
tributions other than the singular isothermal sphere (given 
that the mass ratio of the subject body to the interior mass 

2 

of the host ^ 1). For point like particles, A ^ ^ 1, 
while for extended bodies the Coulomb logarithm is obtained 
from N-body simulat i ons, that give values o f /nA ^2 — 7 
(jSpinnato et al.ll2003l : [Penarrubia et al.ll2004l ). The calcula- 
tion of tfric requires the typical stellar velocity, while we 
are actually measuring the velocity fields of the ionised gas. 
The circular velocity of the stars in the central region of 
disk galaxies will be 1.5-2 times slower than the ionized gas 
(Ganda et al. 2006) due to asymmetric drift. We scale our 
measured velocity values to take that into account. 

We calculate the dynamical friction timescale for clus- 
ters (taken to be point-like test particles) of masses between 
10^ and 10^ M©, that start at an initial distance of ~ 500pc 
to the center. We find that typically clusters with masses 
above 2 x 10^ M© will make it to the dynamical center 
within ^ 2 X 10^ years (see Figure [T]). These findings are 
i n ver y good agreement with the recent study by Bekki] 
(|2010 h. who performed numerical simulations on dynami- 
cal evolution of disk galaxies and investigated the orbital 
evolution of star cluste rs influenced by dynamical friction 
against disk field stars. llBekkil (|2010l ) finds that dynamical 
friction of star clusters against disk field stars is much more 
effective in orbital decay of star clusters in comparison with 
that against galactic halos in disk galaxies. Moreover, dy- 
namical friction seems to be most effective in disks with 
disk masses lower t han 10^ o wing to smaller stellar 
velocity dispersions. 'Milosavljevi^ (|2004h argued that dy- 
namical friction timescales are too long to bring massive 
clusters to the centers of late-type spiral galaxies. However, 
recently lAgarwal &; Milosavlievid (|2010h pointed out that 
iMilosavlievid (|2004h only considered migration from a dis- 
tant location in the disk, and indeed clusters that form close 
to the galactic center can reach the center and merge with 
the nuclear cluster. 



4 DISCUSSION 

The main goal of this paper is to verify whether the NC 
and the KC in late- type, bulgeless galaxies always coincide 
and whether the presence or mass of a NC is related to the 
kinematic state of its host galaxy. As discussed in Section O 
we indeed find that for galaxies with high fidelity velocity 
fields, the NC, the KC and the P C coincide within the errors. 
A similar result was obtained bv i Trachternach et al.l (|2008h 
who conclude, based on HI rotation curves, that there are 
no systematic offsets between optical an d kinematic centers. 
However jTrachternach et aD (|2008h also conclude that. 
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Figure 5. The derived rotation curves for the galaxies in our sample. The lines overplotted to the data points give the best fit to Eq.[T] 
The abbreviation next to the galaxy name indicates whether the galaxy was observed with SparsePak(SP) or with Ppak. 



based on HI observations, non-circular motions are small, es- 
pecially in late Hubble types. As has become evident, this is 
not the case for velocity fields with higher spatial resolution 
and derived from the ionized gas. 

Our derived velocity fields compare well with literature 
results. For example, Ganda et al. (2006) have published H/3 
velocity fields for NGC3346 and NGC3423. While they did 
not perform kinematical modeling of their data, i.e. did not 
publish rotation curves or the KG location, the general shape 
and orientation of their velocity fields agrees well with those 
published here. 

Two of our galaxies, NGG5964 and NGG6509, are 
also part of t he re cent HI study of bulgeless spirals by 
I Wat son et al.l (|2010l V Gomparing our Ra velocity fields to 
their HI data, we find that the position angles and velocity 
scales are in good agreement between atomic and ionised 
hydrogen. The advantage of the HI maps is their large field 



of view, which covers the entire galaxy and thus the com- 
plete rotation curve. Our Ha velocity fields cover a good 
fraction of the rotation curves, but for some galaxies we are 
restricted to the central part of the velocity gradient and do 
miss the turnover of the velocity curve, like in the case of 
NGG 5964. This results in large error bars on the location 
of the KG. We com pare the KG position derived from HI 
(| Wat son et al.l lioiol ) and Ha and find good agreement in 
the case of NGG 6509 but for NGG 5964 the oflFset is larger 
than the measurement error. This is reflected in the fact 
that the best flt to the velocity fleld of NGG 5964 appears 
unreliable. 

We have used our measurements to derive realistic dy- 
namical friction timescales for all sample galaxies (see Sec- 
tion [33]). We flnd that the threshold mass for a star cluster 
to be able to migrate to the center efficiently from a real- 
istic distance (^ 500pc) is 1 — 2 x 10^ M0(in very good 
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Figure 6. Left Projected position of all nuclear clusters, relative to the kinematic centers of their respective host galaxy. Size of the 
crosses denotes the Icr errorbar in the kinematic centers. Right. Offset of the photometric center relative to the kinematic center of the 
non-nucleated galaxies. Here, the size of the crosses denotes the combined Icr error in the kinematic plus photometric centers. The colors 
of the symbols indicate how well the velocity fields can be modelled. Black symbols indicate a model with high fidelity, blue symbols 
denote models with low trust, and red symbols denote galaxies that can not be well modelled due to interactions, or non-circular motions 
(see text for details). 




Figure 7. Dynamical friction time as a function of cluster mass, for a cluster starting 500pc away from the dynamical center of the 
galaxy, and assuming a fixed value of In A = 7. Two example calculations are shown: a galaxy with a nuclear cluster (NGC 4204, left), 
and a non- nucleated galaxy (NGC 5789, right). In both cases, a cluster of ~ 2 x IO^Mq would sink into the dynamical center within two 
Gyrs. 



agreement with the numerical simulations bv lBekkil t01(t )). 
This is inde e d the observed lower mass threshold for NCs. 
iBoker et al] (|2002l ) find that there is a lower luminosity 
cut-off to the NC luminosity function at M/ = —9. For 
a known M/L this w ould imply a lower limit mass cut-off 
for NCs. lWalcher et al. (2005) published dynamically deter- 
mined M/L ratios for a sample of nine NCs with a mean of 



M/L/ = 0.6. This would translate into an observed lower 
mass limit for NCs of ^ 10^ M©. 

We thus find that our data are entirely consistent with a 
NC formation scenario in which a massive seed cluster forms 
within ^ 500pc from the center of the galaxy, and spirals in 
to the center due to dynamical friction. This scenario has the 
advantage of automatically accounting for galaxies without 
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a nuclear star cluster, implying merely that no suitable seed 
clusters formed close enough to the center in these disks. 
We also note that this lower mass threshold is one order 
of magnitude lower than the threshold mass of - 10^ Mo 
derived by El amm-Altenburg Kroupal (|2009l ) for efficient 
accretion of gas into the NC, which would lead to repetitive 
bursts of star formation there. Therefore, it may be that 
further NC growth through in-situ star formation starts out 
quite slowly, but picks up in pace as the NC grows. Mass 
growth may be further supported by accretion of other mas- 
sive inspiralling clusters. Soon after its formation the seed 
cluster would then satisfy our definition and observations of 
a nuclear star cluster, namely being massive, sitting close 
to the center, having undergone recurrent star formation, 
and being compact. This formation and growth mechanism 
is also supported by the fact that NCs are observed to rotate 
(|Seth et al.ll2QQ8l . l2010l ). 

The long dynamical friction timescales we infer for mas- 
sive star clusters that are initially located at significant (> 1 
kpc) distances from the galaxy center do not favour mod- 
els where the NC grows only through accretion of globular 
clusters. We emphasize on the other hand that we have not 
found a confirmed case where the NC sits outside of the KC. 
We thus cannot rule out in-situ formation scenarios for NCs 
from our data. 



5 CONCLUSIONS 

We have presented two-dimensional Ha velocity fields for 20 
late-type, disk-dominated spiral galaxies, the largest sample 
of bulgeless disks with high-resolution Ha velocity fields to 
date. 

We fitted the data with kinematic models in order to 
derive rotation curves and the location of the dynamical cen- 
te rs. Most r o tation curves are well-fit by the arctan form 
of ICourteaul (|l997l l used for earlier- type spirals. We find 
that the velocity fields span a broad range of morpholo- 
gies. Some galaxies show regular rotation, which allows ac- 
curate determination of the kinematic center. However, only 
2 out of 20 galaxies have a completely regular velocity field. 
Many galaxies without bulges, but with strong bars show 
steep rises of their "velocity gradients" ("rotation curves") 
in the inner part. These steep rises are not necessarily due 
to mass concentrations, but to streaming motions along the 
bar. However, quite a few have some degree of irregular gas 
motion, which in nearly all cases can be either attributed 
to the presence of a bar or is connected to a rather patchy 
appearance of the Ha emission and the stellar light. Thus 
most galaxies in the sample show strong gas motions that 
cannot be attributed to the overall gravitational potential 
of the galaxy, implying l ow surface mass densities (compare 
IPalcanton Stilpl 20 id) and implying that many bulgeless 
galaxies are not in dynamical equilibrium. 

There appears to be no systematic difference in the kine- 
matics of nucleated and non- nucleated disks. For galaxies 
with regular, well-sampled velocity fields, the photometric 
center, the nuclear cluster and the kinematic center coin- 
cide within the errors. These centers also coincide for quite 
a few of the not-so-regular galaxies (in total for 13 out of our 
sample of 20 galaxies). However, we also find that nuclear 
clusters also occur in galaxies with disordered rotation fields. 



Hence, the large-scale velocity field is not a good predictor 
for the presence or mass of a nuclear cluster. Many formation 
scenarios for nuclear clusters invoke off-center cluster forma- 
tion and subsequent "sinking" of clusters due to dynamical 
friction. We confirm that this scenario is viable for clusters 
that form within 500pc of the center of the galaxy, as 
the dynamical friction timescales inferred from our data are 
consistent with this scenario. More distant globular clusters 
do not seem to be promising candidates for nuclear cluster 
seeds, due to their long dynamical friction timescales. On 
the other hand, we point out that we cannot rule out an 
in-situ formation scenario. 
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